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Abstract 

Non-Hermitian Hamiltonians H ^ possess the real (i.e., observable) spectra 
inside certain specific, "physical" domains of parameters V = V{H). In general, 
the determination of their "observability- horizon" boundaries dT> is difficult. We list 
the pseudo-Hermitian real by matrix Hamiltonians for which the "prototype" 
horizons dV are defined by closed analytic formulae. 



1 Introduction 



In Landau's textbook [Ij on Quantum Mechanics the emergence of an instabihty in 
a system is illustrated via a particle in the potential V{x) = G/\x\'^. The critical 
value G{rnin) = — 1/4 of its strength represents a "horizon" beyond which the particle 
starts falling on the center. Vice versa, the system remains stable and physical on all 
the interval V = (—1/4, oo) of couplings G. From the pragmatic point of view the 
Landau's example is not too well selected since the falling particle should release, 
hypothetically, an infinite amount of energy during its fall. A better example of the 
loss of stability is provided by the Dirac's electron which moves in a superstrong 
Coulomb potential: In the language of physics, particle-antiparticle pairs are created 
in the system beyond a critical charge (^(max) = 137 in suitable units |2]). 

What is shared by the above two sample Hamiltonians H{X) is that they are well 
defined in a certain domain V of parameters A while they lose sense and applicabil- 
ity for parameter (s) beyond certain horizon(s) X{max/min)- On a less intuitive level, 
similar situations have been studied by Kato |3j. He considered certain finite- matrix 
toy Hamiltonians H{X) in complex plane of A and deduced that the related (in gen- 
eral, complex) spectra -En (A) change smoothly with the variation of the parameter 
A unless one encounters certain critical, "exceptional" points A^^^-*. 

Several recent microwave measurements [1] confirmed the observability of the ab- 
stract Kato's exceptional points A^^^-* in practice. These experiments re-attracted 
attention to the theoretical analyses of the EP horizons, say, in nuclear physics where 
many nuclei can, abruptly, lose their stability [5, 6j. The presence of EPs may be 
also detected in the random-matrix ensembles with various interpretations [7J and in 
optical systems (where EPs are called "degeneracies" [8j). In classical magnetohy- 
drodynamics the Kato's exceptional points may even happen to lie inside the domain 
of acceptable parameters, separating merely the different dynamical regimes of the 
so called a^— dynamos [9J. 

The obvious theoretical appeal of the problem of stability may be perceived as 
one of the explanations of the recent growth of popularity of the so called pseudo- 
Hermitian Hamiltonians in quantum physics [10]. Indeed, one of their distinctive 
features is that their spectra are real (i.e., observable) in parametric domains V 
with, sometimes, very complicated and strongly Hamiltonian-dependent shape of 
their EP boundaries &D. For an uninterrupted development of their study it was 
very fortunate that a successful semiquantitative description of the spiked-shaped 
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horizons dV has been found, by Dorey, Duncan and Tateo [TT], for an important 
class of analytic (polynomial and power-law, often called PT— symmetric) potentials 
A) with promising relevance in quantum field theory [121 fT3] - 
In a few of our own recent studies of EPs in pseudo-Hermitian Hamiltonians 
H [H] - [21] we paid detailed attention to the possibilities of a deeper geometric 
understanding of the structure of the domains V (H) of quasi-Hermiticity (the name 
means that the spectrum remains real for parameters inside T> - cf. ref. |5] for an 
older, well written introduction of this concept). After we review some of the known 
results in section [2] we shall combine, in section [3l the methods of algebra (of solvable 
equations) and analysis (of elementary curves) in a new approach to the problem. In 
this way, the list of results of ref. [18] (based mainly on computer-assisted symbolic 
manipulations) and of ref. [19] (which used, predominantly, pertubation-expansion 
methods) will be complemented by a number of new non-perturbative items. They 
will be discussed and summarized in sections H] and O 

2 Matrix models 

2.1 Inspiration: two-dimensional Hilbert space 

The first nontrivial schematic illustration of the current Schrodinger's bound-state 
problem is provided by the two-by-two real-matrix model 



In the context of PT— symmetric Quantum Mechanics [T2l [22] , the parallel two-by- 
two example is very similar 



It can still be considered Hermitian (or, in the language of the review paper [5], 
quasi-Hermitian) after one redefines the Hilbert space accordingly (cf. [H] for more 
details). 




Its three-parametric spectrum is always real and, therefore, observable. 
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In the language of phenomenology, one notices an important complementarity 
between the parameter-dependence of the two toy spectra E±{a, b, d) and E'^{a, b, d). 
In the "classical", former example, all of the energies E±{a,b,d) remain safely real. 
The second, primed model is less easy to deal with. There exists the whole set of 
the eligible two-by-two metric operators = 0^ > which define the inner product 
in the corresponding two-dimensional toy Hilbert space Ti' (cf. [U]). Thus, in spite 
of its manifest non-Hermiticity in the auxiliary two-dimensional Hilbert space Ti. 
(where the metric is the Dirac's simplest identity operator), the operator H'{a,b,d) 
represents an observable and remains safely compatible with all the postulates of 
Quantum Mechanics (cf. reviews O [12] for more details). 

For H = H{a, b, d) = the three-dimensional physical domain I^(a, b, d) of 
parameters giving real spectra coincides with all M^. In contrast, for each individual 
choice of the parameters a, b and d, the quasi-Hermiticity property of the primed 
Hamiltonian H'{a, b, d) must be guaranteed and proved. In general, the reality of the 
bound-state energies E'^{a,b,d) and/or the stability of the primed system can only 
be achieved inside a perceivably smaller domain T''(a, b, d) with the easily specified 
EP horizon. 



Thus, the interior of the non-compact manifold T>'{a, b, d) is specified by the single 
elementary constraint b G (— |a — d\, \a — d\). 

Of course, for qualitative considerations the variability of parameters a and d is 
entirely redundant. It makes sense to get rid of them by the multiplicative re-scaling 
of all the parameters and by the subsequent shift of the energy scale leading to the 
"generic" choice of a = — 1 and d = 1. An extension of this argument has been 
formulated in refs. [16j - [19j. A support has been found there for the study of the 
very special family of matrix models H^^\ In our present paper we just intend to 
add new results showing that and how the respective quasi-Hermiticity domains V 
can be described, at the lowest dimensions, by non-numerical means. 

2.2 Tridiagonal chain models 

One of the most surprising byproducts of our studies [H] - [21] was an empirical ob- 
servation that certain particularly simple (i.e., in the sense of regularity, "canonical") 
shapes of the J— parametric quasi-Hermiticity domains T'(^) can be found for a class 
of the generic, maximally simplified pseudo-Hermitian Hamiltonians H^^^ chosen in 
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the following A^— dimensional and tridiagonal "self-dual" 
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with a J— plet of real couplings A = {gi, (72, ... , gj) and with the dimensions N = 2J 
or N = 2J + 1. In a more explicit formulation, at any dimension we found the 
coordinates of all the maximal-coupling spikes of the horizon dV^'^^ in closed form, 
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{spike) 



±{N -n)n, 



n 



1,2, 



,J. 



(2) 



Although this result looks easy, its derivation from the underlying algebraic equations 
required extensive computer-assisted symbolic manipulations and nontrivial extrap- 
olation guesswork |18j . Moreover, this closed-form description of the positions of the 
protruded spikes of the horizon dV^^'' (called, in [TB], "extremely exceptional" points, 
EEPs) has been complemented, in our subsequent paper [19], by the strong-coupling 
description of dV^^^ based on perturbation ansatz 



9n 



7„(t) = t + t2 + 



(3) 



This formula extrapolated, to all J, the rigorous J < 2 fine-tuning rules as derived 
in refs. [TB]. A posteriori, using sufficiently small "redundant" parameters t <^1, 
it proved valid in nonempty open vicinities of all the EEP vertices. 

At the larger ts, i.e., far from the EEP spikes, the determination of the physical 
horizons dV^^^ of our models H^^^ with J free real parameters becomes a more or 
less purely numerical task at the higher Js [20] ■ Up to now, non- numerical exceptions 
with N = 2 and A^ = 3 have been reported in [T3] (where the very easy localization 
of the one-parametric interval P*^^^ = (—1, 1) of (71 has been made) and in [18] 
(mentioning the very similar result P*^^^ = (— -\/2, -\/2)). For the next two dimensions 
A^ = 4 and A^ = 5 with two parameters, the explicit construction of the planar 
curves can be also found in ref. [18|. In what follows we intend to complement 

and extend these observations beyond J = 2 and to show that the closed-form 
constructions of the prototype horizons dV^^^ remain feasible up to the dimension 
as high as A^ = 11. 
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2.3 Secular equations 

Once we choose = 2 J or N = 2J + 1, abbreviate = s and, at all the odd 
dimensions = 2J + 1, ignore the persistent and trivial "middle" energy level 
j^{2j+i) _ ^-j^gj ^YioX all the secular equations det (^H^^^ — E^ = have 

the same polynomial form, 

j .^-p+(|^ j .^-g-(|^ j s^-^i? + ... = o. (4) 

At all J and A^, the coefficients P,Q,R,... are real polynomial functions of the 
squares gl of the J— plets of our real matrix elements. Once all the energies are 
assumed real (i.e., equivalently, once all the roots Sk of eq. (jl]) happen to be non- 
negative), we immediately deduce the following relations tractable as necessary con- 
ditions imposed upon our coefficients in (jlj). 




In the opposite direction, the set of the necessary inequalities P > 0, Q > 0, ... is in- 
complete as it does not provide the desirable sufficient condition of quasi- Hermiticity. 
It admits complex roots s in general (take a sample secular polynomial (s^ + l)(s — 2) 
for illustration). 



3 Domains P^^J) ^nd p(2J+i) 

Obviously, for a given prototype Hamiltonian H^^^ and under the constraints ([5]), 
the determination of the quasi-Hermiticity domain T'(^) = V (^H^^^^ is equivalent 
to the guarantee of the non- negativity of all the J roots Sk of eq. (jlj). The explicit 
forms of the corresponding sufficient conditions will now be given for the ffist ten 
smallest matrix dimensions A^ = 2,3,...,ll. 
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3.1 Non-negativity of the root of eq. (jll) at J = 1 

At J = 1 the linear version s — P = of secular eq. (jlj) has the single root Sq = P. 
The non-negativity of this root is equivalent to the non-negativity of the coefficient 
P. This means that in terms of the single coupling gi = a available at J = 1, the 
necessary and sufficient criteria of the observability of iJ*-^-* = H^'^\a) or 
H^^\a) read P^'^\a) = 1 - > and P^^\a) = 4 - > Q, respectively. Thus, 
in a way transferable, mutatis mutandis, to any dimension, the explicit definitions 
V^'^\a) = (—1,1) and T)^^\a) = (— -\/2, v^) of the quasi-Hermiticity domains may 
be re-read as definitions of the corresponding EP horizons cfD^'^\a) = {—1, 1} and 
dV(^){a) = {-V2, 

3.2 Non-negativity of all the roots of eq. ([4l) at J = 2 

At J = 2 the quadratic version s^ — 2Ps + Q = 0of secular eq. (jlj) has two roots 
s± = P ± \/P'^ — Q. These two roots remain real if and only if i? = — Q > 0. 
In the subdomain of parameters where P > they remain both non-negative if and 
only if P > and Q > 0. We can summarize that the required sufficient criterion 
reads 

p > , p^ > g > . (6) 

In an alternative approach, without an explicit reference to the available formula for 
s±, let us contemplate the parabolic curve y{s) = s"^ — 2 P s which remains safely 
positive, in the light of our assumption ([5]), at all the negative s < 0. This means 
that this curve can only intersect the horizontal line z{s) = —Q at some non-negative 
points s > 0. 

In this way the proof of non-negativity of all the roots of our secular equation 
degenerates to the proof that there exist two real points of intersection of the J = 2 
parabola y{s) with the horizontal line z{s) (which lies below zero) at some s > 0. 
Towards this end we consider the minimum of the curve y{s) which lies at the point 
So such that y'{sQ) = 0, i.e., at sq = P. This minimum must lie below (or, at worst, 
at) the horizontal line of z{s) = —Q < 0. But the minimum value of y{so) is known, 
y{P) = — P^. Thus, the condition of intersection y{so) < z{so) gives the formula 
P2 > Q. QED. 

Marginally, it is amusing to notice that once eq. (jS]) holds, the inequality P^ — Q> 
is equivalent to the reality of the roots simply because P^ — Q = [si — S2)^/4. In 
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fact, even for some other two-parametric matrices, precisely this type of requirement 
is responsible for an important part of the EP boundary dT> (cf. refs. [T7] for 
details). 



3.3 Non-negativity of all the roots of eq. ([4l) at J = 3 

Neither at = 6 nor at = 7 the sufficient condition of non-negativity of all the 
energy roots s is provided by the three necessary rules P > 0, Q > and -R > 
of eq. Let us return, therefore, to the second method used in paragraph 13.21 
and derive another, "missing" inequality needed as a guarantee of the reality of the 
energies. In the first step one notices again that all the three components of the 
polynomial 

y{s) = -3Ps^ + 3Qs = R, J = 3 

remain safely non-positive at s < 0. Whenever the roots are guaranteed real, their 
non-negativity s„ > with n = 1,2,3 is already a consequence of the three con- 
straints ([5]). The necessary condition of their reality is less trivial but it still can be 
deduced from the shape of the function y{s) on the half-axis s > 0, i.e., from the 
existence and properties of a real maximum of y{s) (at s = and of its subsequent 
minimum (at s = At both these points the derivative y'{s) = 3s^ — 6Ps + 3Q 
vanishes so that both the roots s± = P ± \JP^ — Q of y\s) must be real and non- 
negative. This condition is always satisfied for the real roots of y(s) since 



5 = p2 _ g ^ 1 [(^^ + S2 - 2 S-^f + (S2 + S3 - 2 S^f + (^3 + Si - 2 S^f 

54 L 



> 0. 



In the next step, the necessary guarantee of the reality of the roots will be under- 
stood again as equivalent to the doublet of the inequahties y{sJ) > R and ?/(s+) < R. 
Here we may insert = 2Ps± — Q and get the two inequalities which are more ex- 
plicit, 

2(p2_g)s_ < PQ-R < 2{P'^-Q)s+. (7) 
They restrict the range of a new symmetric function of the roots, 

PQ-R = ^ [S1S2 (si + S2 - 2 S3) + S2S3 (s2 + S3 - 2 si) + S3S1 (S3 + si - 2 S2)] . 

After another insertion of the known s± we arrive at a particularly compact formula 

2(p2 - g)3/2 > R-3PQ + 2P^ > -2 (P2 - Q)3/2 
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or, equivalently, 

4 _q)^ > (^R~3PQ + 2P^y . 

Due to the numerous cancellations the latter relation further degenerates to the most 
compact missing necessary condition 

3P^Q^ + 6RPQ>4Q^ + R^ + 4RP^. (8) 

Our task is completed. In combination with eqs. ([5]), equation ([8]) plays the role of 
the guarantee of the reality of the energy spectrum. 

3.4 Non-negativity of all the roots of eq. ([4l) at J = 4 

In a search for the non-negative roots of the quartic secular equation 

det -EJ) =x^ -APx^ + QQx'^ -ARx + S = y{x) + S = (9) 

we note that all the four A^— dependent coefficients P, Q, R and S again evaluate as 
certain polynomials in the squares of the four coupling parameters Qk, k = 1,2, 3, 4. 
Once all these four expressions are kept non-negative, the curves y{x) and z{x) = —S 
do not intersect at x < 0. At a; > they do intersect four times at a; > (as required), 
provided only that the three extremes of y{x) can be found at the three non-negative 
real roots xi^2,3 of the extremes-determining equation 

y\xi,2,3) = 4 (a;?,2,3 - 3 P a:?,2,3 + 3 Q x,,2,3 - R) = . (10) 

In an ordering < Xi < X2 < X3 of these roots we arrive at the three sufficient 
conditions 

y{x^)<-S, y{x2)>-S, y{x3)<-S (11) 

guaranteeing that the parameters lie inside V^^^ or T>^^\ 

All the three quantities x^ satisfy the cubic equation y'{x) = so that its pre- 
multiplication by x enables us to eliminate the three fourth powers from y{x), 

^1,2,3 = 3 P a:^ 2,3 — 3 Q 2,3 + ^1,2,3 • 

Their insertion reduces all the three items of eq. (ITT]) to the other three intermediate 
polynomial inequalities of the third degree, 

-P X? 3 + 3 Q x^ 3 - 3 i? xi,3 + S <0, 



9 



-P xl + 3Q xl - 3 Rx2 + S > 0. 

Repeating the same trick once more, the ehmination ofxfgs = 3Px^23~3Q Xi^2,3 + 
R gives an equivalent triplet of inequalities 

-Bxl + 2B^^^Cxi<B^D, (12) 
-Bxl + 2B^/^Cx2>B^D, (13) 
-Bxl + 2B^/^Cx3<B^D. (14) 

Here, the old abbreviations B = P"^ - Q and 2 B^/'^ C = PQ - R plus a new one, 
3 B^ D = PR — S enable us to re-scale xi^2,3, = y/BYi,2,3 which yields our final 
triplet of quadratic-equation conditions 

-2CYi + D >0, (15) 
Yi - 2CY2 + D <0, (16) 
Y^ -2CY3 + D>0. (17) 

The auxiliary roots Y± = C ± \JC'^ — D must be real and non-negative. In this way 
we must guarantee that D > and > D. The conclusion is that eqs. ( |T5i) - ( |T71) 
degenerate to the four final elementary requirements 

Yi<Y^<Y2<Y+<Y;. (18) 

They complement the inequalities i? > 0, Q > and — 1 < C — y^l + QZ-B < 1 and 
form the complete algebraic definition of the domains "D'-^-* and V^^\ 

3.5 Non-negativity of all the roots of eq. ([4l) at J = 5 

Let us finally proceed to H^^^ with = 10 and/or = 11 which leads to the 
"unsolvable" secular equations of the fifth degree, 

x^ - h P x^ + IQQ x^ - Rx^ + h S X - T = y{x)-T = 0. (19) 

From our present point of view the problem of the construction of the respective 
horizons dV^^^ remains solvable exactly since the derivative y'{x) is still of the mere 
fourth degree, 

-y'{x)=x^-4:Px^ + 6Qx'^-4:Rx + S. (20) 
5 

The exact, real and non- negative values xi < X2 < X3 < X4 of the four roots of y'{x) 
may still be considered available in closed form. 
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In a way which parallels our preceding considerations we may assume that the 
five A^— dependent non-negative coefficients P>0,Q>0, R>0, S>0 and T > 
obey also all the additional inequalities derived in the preceding sections. In a more 
detailed description, we may then treat our secular problem ( IT9|) as a search for the 
graphical intersections between the (nonnegative) constant curve z{x) = T and the 
graph of the polynomial y{x) of the fifth degree (which can only be nonnegative 
at X > 0). Inside the domain T)^^\ the quintuplet of the (unknown but real and 
nonnegative) physical energy roots Xa, Xb, Xc, Xd and Xe may be assumed compatible 
with the obvious intertwining rule 

^ ^ Xa < Xi < Xh < X2 < Xc < X^ < Xd < < Xe ■ 

The way towards the sufficient condition of the existence of the real energy spectrum 
remains the same as above, requiring 

l/(xi)>T, y{x2)<T, yixs)>T, y{x^) <T . (21) 

The lowering of the degree should again reduce eq. (pTj) to the quadruplet 

w(ri)<0, w(F2)>0, w(l3)<0, w(F4)>0. (22) 

where the re-scaling 0:1^2,3,4 = ^1,2,3,4 applies to the arguments of the brand new 
auxiliary polynomial function of the third degree in F, 

w{Y) = Y'^ -^CY"^ + ^DY -G . 

Besides the same abbreviations as above, we introduced here a new one, for PS—T = 
4 5^/^ G. The new and specific problem now arises in connection with the necessity 
of finding the three auxiliary and, of course, real and non-negative roots of the 
cubic polynomial w{Y). Once we mark them, in the ascending order, by the Greek- 
alphabet subscripts, we should either postulate our (in principle, explicit) knowledge 
of their real and nonnegative values Y^ <Yi3 <Y^ or, in another perspective, we have 
to add the above-studied conditions which restrict the range of the three coefficients 
G, D and G in the cubic polynomial w{Y). 

This being said, the rest of the story is easy to tell. Once we parallel the same ge- 
ometric argument as used in our previous subsections, we may immediately conclude 
that our "last feasible" specification of the domains T)^"^^^ and T)^^^^ will be given by 
the following set of the inequalities, 

Yi<Y,,<Y2<Yp<Y:,<Y^<Y^. (23) 
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This is the last algebraic formula which defines the domains and T)^^^\ Any 

extension of the recipe beyond = 11 would suffer from the necessity of using 
certain purely numerically defined auxiliary functions of couplings Qk- 



4 Discussion 

4.1 Reparametrizing the couplings 

The existence of the algebraic formulae which determine all the boundaries &V^^^ 
up to = 11 opens a way towards a non-perturbative extension of the results of 
refs. [18] and [19] beyond the strong-coupling dynamical regime. In such a setting, 
the old perturbation ansatz ([3]) can be re-interpreted as a precise, non-perturbative 
change of variables — ^ Gk- During such a process, the redundant value of t 
may be fixed arbitrarily (cf. the construction of the planar curve dT>^'^'> in section 
3.1 of paper [19] for illustration). Due to the refiection symmetries of our hedge- 
hog-shaped horizons d'D^^\ one is also allowed, without any loss of generality, to 
restrict attention to the subdomain of T'(^) with positive gkS. In such a setting the 
new, rescaled real couplings 7j^'' = a, "fj^i = (3, . . . should remain non-negative 
and smaller than one, 'y^^^ G (0,1). For illustration one may recollect the most 
elementary two-by-two Hamiltonian 



with the two-point spectrum E^' = ±^/a . With respect to the new parameter we 
have to set V^^\a) = (0,1) since there are no additional constraints. 

4.2 New forms of approximations 

We should emphasize that in comparison with our previous work, the most signif- 
icant progress has been achieved here in the new non-approximate and complete 
description of the structure of the boundaries of the domains T)^^^ at N = 6 and 
A^ = 7 by means of our key inequality (I7j). In order to stress some of its merits, let 
us now add a few comments on this form of the rigorous guarantee of the reality of 
the energies at J = 3. 

Firstly, let us set = B + Q and switch just to the postulates B > and Q > 0, 
re- classifying the expression P = +y/B + Q itself as a mere formal abbreviation. 




a e (0,1) 



(24) 
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This enables us to insert 



= ^5 + Q ± > 
in eq. ^ which prescribes s_ < C < s+ with PQ - R = 2 B'^/^ C > 0, i.e., 

Q 



l + q-l < C < Jl + q + l, g=^G(0,oo). (25) 

Once we notice that PQ = Q \JB + Q = qB^I'^^X + g we may return from the 
auxihary C to the original R and rewrite eq. (!25|) as a perceivably simplified one- 
parametric constraint 

imposed upon the rescaled polynomial R. It completes our specification of the phys- 
ical domain D^^'''') in terms of the coefficients -B > 0, Q > and R. We see that 
with 

1 + _ /TT^ = 3 2 _ 1 _9_ _ ^ . 
V2 /V 8^ 128^ 64^ 1024^ 

the two-sided inequality (l26l) is fairly restrictive in the i?— direction, especially at the 
smallest ratios q = Q/B. 

4.3 Pair wise confluences of the levels 

Several aspects of the "first nontrivial" J = 3 problem have been skipped in the 
main text but they definitely deserve more attention in the discussion. For the sake 
of definiteness let us choose just = 6 and use 



9i = c= ^5{l-j), g2 = b = 2^2{l-P), g, = a = 3VT^ 

with parameters a,/3,7 G (0, 1) entering the six-by-six matrix ([1]). Its secular equa- 
tion (jlj) determines the spectrum which remains real, by definition, whenever all the 
couplings lie inside the domain T)^^\ It is easy to deduce that the latter domain is 
circumscribed by the ellipsoidal surface given by the equation 

P = - (a2 ^ 2 6^ + 2 - 35) /3 = . 

The other two obvious constraints read 

3 g = 6^ + 2 c^a^ -Ub'^ + 28 - 34 + + 259 + 2 b'^c^ > 
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and 

-R = a^c^ - 10 b'^c^ + 30 c^a^ + 225 - 30 - - 25 6^ - 225 - 150 6^ > . 

The last constraint needed to define V^^^ is then given by eq. ([7j). In its hght one can 
spot certain new structures in dV^^^ where, for example, the total EEP confluence 
of all the energies can be preceded by some incomplete, pairwise coincidences among 
the six levels in question. For example, an "innermost" pair of the energies can 
coincide at = while, independently, the other two "outer" doublets of the energies 
are allowed to coincide at the two symmetric non-vanishing values E = ±4z at 
an unknown parameter z G (0,1). Alternatively, the two "outermost" levels [with 
s = Smax = 5?/ where y G (0, 1)] can stay real while the confluence only involves the 
two internal energy doublets at a shared value of s = Smin = 4x where x G (0, 1). 
In the latter scenario one has to reproduce the two-parametric relation 

(s-s^ax) [s- s^inf = s^- (32 + 25 y^) s^ + (256x^ + 800a;^?/2) s-6400xV = 

which, implicitly, deflnes the third sub-surface. Let us concentrate on the former, 
slightly simpler scenario where the surface of the three-dimensional domain T)^^^ can 
be visualized as composed, locally, of the two eligible smooth sub-surfaces which 
intersect along a certain "double exceptional point" (DEP) curve. 



4.4 Pairwise confluences of exceptional points 



In terms of the single free parameter z of the latter particular scenario, the DEP secu- 
lar equation degenerates to the formula E"^ {E + 4z)^ {E — Az)^ = 0, to be obtained, 
in the one-parametric DEP limit, from eq. (jlj), i.e.. 



s - (4z) 



s-^ -2{Azys' + {Az^s = 0. 



(27) 



It is fortunate that the necessary analysis can still be performed non-numerically since 
equation (1271) is easy to compare with the true secular equation (jl]) with coefficients 
given in paragraph 14. 3[ As long as the factorizable coefficient at s° must vanish, we 
get the first DEP constraint 



ac 



-M5a-F (l5 + 5 6^ 
so that we may eliminate 



ac^ 15 a - (l5 + 5 6^ 

l5 + c^ + 5b'^ 
c2 + 15 
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In the quadrant oi a — b — c space with positive a the plus sign must be chosen, 

" = ' + ^^ 
i.e., we have 3 > a > 1 in the closed formula for 

b^ = l (c^ + 15) (a - 1) (28) 
5 

or, alternatively, for 

c = 15 . 

a — 1 

This result is to be complemented by the other two relations 

3Q{c,b,a) = 32z'^ , R{c,b, a) = 128 z\ 

A straightforward elimination of gives the second DEP condition 

-66 - 36 6^ + 4 c^a^ - 189 + 252 0^-4 b^a^ - = 

with the two compact roots 

a| = 2 - 33 - 2 6^ ± 2 Vc^ + 30 - 2 b^c^ + 225 + 24 62 + . 

The acceptable one must be non-negative. For this would mean that 2c2 > 
33 + 2 6^ while, at the same time, 63 + 12 6^ > 84 c^. These two conditions are 
manifestly incompatible so that we must accept the upper-sign root a\ which is 
automatically positive for all the large 2 > 33 + 2 6^ and which remains positive 
for all the smaller constrained by the requirement 

84 > 63 + 12 6^ . 

After the insertion of the definition (!28|) of 6^ we arrive at the condition 

12 

84 > 63 + y (a - 1) (l5 + c^) (29) 
with a non-empty domain of validity. 

5 Conclusions 

According to the abstract principles of Quantum Mechanics, observable quantities 
(say, the energies Eq < Ei < . . . of bound states) should be constructed as eigen- 
values of certain self-adjoint operators of observables (i.e., in our illustration, of a 
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Hamiltonian H acting in some physical Hilbert space of states Ti). Of course, an 
explicit representation of Ti can prove complicated. Hence, the idea emerged of an 
introduction of a simpler, "auxiliary" Hilbert space (to be denoted here as Ti*^"""^^). 

Although the origin of the latter idea can be traced back to the very early days 
of Quantum Theory [13], the feasibility of its separate implementations have long 
been treated as a mere mathematical curiosity (cf., e.g., [23] for illustration). Among 
the people who felt forced to take it really very seriously were the nuclear-physics 
specialists. Incidentally, in their so called IBM models of atomic nuclei, one of the 
"natural" (often called Dyson's) fermion-to-boson mappings Ti ^ q-i^aux) happened 
to simplify the solution of the Schrodinger equation considerably. In 1992 the topic 
has nicely been reviewed by Scholtz et al [5]. 

In 1998, many other physicists of different professional orientations (ranging from 
supersymmetry ^25] and field theory |26j to cosmology [27] etc) have got involved 
when Bender and Boettcher |28j attracted their attention to the specific class of 
examples 

(x) + l/^'*"^) (x) ?/'(''"^) (x) = E V-^""^) (x) (30) 

2m dx^ 

(cf. also refs. [221 ED]) where the potentials may be complex but where the spectra 
remain real [TT]. Thus, although the physical space H, itself proves complicated 
(mainly, due to its highly nontrivial definition of the inner product [31]), the auxiliary 
space can often be chosen in its most common representation Ej2{M) of the 

square-integrable complex functions of one variable. In this setting (cf. its recent 
review [12]), the difference between the clarity of the properties of the space of 
functions '}{(°''^^) = Ej2{Si) and the perceivably more complicated character of the 
physical states in Ti emerges due to the nontriviality of the "physical metric" Q ^ I 
in Ti (cf. ref. |f51| for a compact review of the necessary mathematical properties of 
this operator). In the context of eq. (I5U]) (possessing also numerous exactly solvable 
special cases [30j), the metric 6 is assumed constructed as the product of certain 
operators V (= parity) and C (= "charge"). 

Beyond the specific class of the differential phenomenological models fl5U]l . another 
obvious choice of Ti^""^) would be finite-dimensional [SIE2]. One of the key merits of 
such an alternative extension of the class of the "tractable" models of bound states 
lies in the fact that the key mathematical proofs of the reality of their spectra may 
become decisively simpler. It is in the latter context, with dim Ti*^""^-' = N < oo, 
where our attention has been attracted to the various specific models where one can 
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extract more information about the shape of the spectral-reahty domain V of variable 
parameters in the Hamihonians. 

As we emphasized in [20] , a nontrivial relationship may exist between the matrix, 
finite-dimensional models and the so called quantum catastrophes interpreted as 
changes of some parameters Aj, j = 1,2,..., J in the Hamiltonians H{X) which 
lead to the loss of the reality (i.e., of the observability) of certain energies En^{X) in 
the spectrum. In a step towards a typology of such a phenomenon we showed here 
that for our PT— symmetric chain models H^^^ of dimensions = 2,3, . . . , 11 a 
non-numerical knowledge becomes available about the parametric domains p(^) in 
which all the energies remain real. This means that the "admissible" J— plets of the 
coupling constants have been shown defined by certain non-numerical means. In this 
way, the explicit control of the stability of the system becomes mediated via purely 
algebraic constraints imposed upon the controllable parameters. 

In the conclusion we feel encouraged to express our belief that many qualita- 
tive and geometric features of the observability horizons dT> assigned to any given 
PT— symmetric Hamiltonian H may be expected to survive its reduction to a series 
of by approximations of the prototype form H^^\ Such a feature would really 
enhance the relevance of our present study of the peculiar self-dual PT— symmetric 
Hamiltonians if*-^-* ^ (^11'-^^^. In particular, it would also assign a deeper meaning 
to our present detailed analysis of the boundaries dV^^^ at the first few dimensions 
up to A^ = 11. In such a setting and perspective it is also appropriate to re-emphasize 
the two reasons of the relevance of our knowledge of the physical domains V. Firstly, 
their boundaries &D are marking the breakdown of the reality and observability of 
the spectrum {En}- Secondly, all the vicinity of these boundaries also represents a 
region where the matrices H^^^ cease to be diagonalizable. Thus, it is the simulta- 
neous degeneracy of the energies and of the wave functions which gives the full and 
deep physical meaning to this horizon of the dynamical stability of the system. 
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